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ff^i Abstract 



We study the effects of breaking S3 symmetry in the neutrino mass matrix for the 
masses and mixing matrix of neutrinos. At zeroth order the model gives degenerate 
neutrino masses and accommodates tribimaximal mixing. We introduce perturbations 
in terms of a small and complex parameter. The perturbations are introduced in a 
manner such that the S3 symmetry is broken by its elements in the same representa- 
tion. Successive perturbations introduce mass splitting, sizable non-zero reactor mixing 
angle and CP violation. This scheme of breaking 5*3 symmetry can reproduce a rela- 
tively large reactor mixing angle as suggested by the recent T2K results. The effective 
neutrino mass is predicted to be large which is testable in the ongoing and forthcoming 
neutrinoless double beta decay experiments. 



1 Introduction 



After the discovery of neutrino oscillations, there has been considerable progress in determin- 
ing the values of the neutrino mass-squared differences and the mixing angles which relate 
the neutrino mass eigenstates to the flavor eigenstates. Recently, T2K [1J, MINOS [2 J and 
Double Chooz [3] experiments have given hints of a relatively large 1-3 mixing angle. This 
has given a fresh impetus to the construction of neutrino mass models which can accom- 
modate a non-zero value of #13. Recently many papers have appeared which reproduce the 
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relatively large value of the reactor mixing angle [I] . The mixing pattern in the quark sector 
comprises of small mixing angles while in the lepton sector the mixing is described by two 
large angles and one small angle. Discrete symmetries have been extensively used in the 
past to obtain two large and one small mixing angle in the lepton sector. The most studied 
phenomenological Ansatz for the neutrino mixing matrix arising from discrete symmetries 
was proposed by Harrison, Perkins and Scott [5] known as tribimaximal mixing (TBM) is 
given by 
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Non-Abelian flavor symmetries [6] have been invoked to account for TBM. With the recent 
T2K results various neutrino mass models based on TBM need to be suitably modified to 
accommodate a relatively large (9 13 . The smallest discrete non-Abelian group is S3 which 
is the permutation group of three objects. A large number of papers j?l |8] have presented 
detailed models based on S3 symmetry. The permutation matrices in the three dimensional 
reducible representation are 
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where the matrices in each equation belong to the same class of S3 
neutrino mass matrix M u under S3 requires 

[S,M v ] = 



1 1 \ 

1, (4) 
Vo 1 0/ 

The invariance of 
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where S is any of the six permutation matrices given in Eqs. (2, 3, 4). The neutrino mass 
matrix invariant under S 3 is given by 



M„ = aI + bD 
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where a and b are, in general, complex and D is called the democratic matrix. In this work 
we consider deviations from S3 invariance of the neutrino mass matrix. We introduce the 
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perturbations which break this symmetry in such a way that the perturbed neutrino mass 
matrix still satisfies the magic symmetry [9| ITU] . for this we choose the perturbation matrices 
to be the S3 group matrices. Indeed any linear combination of these matrices will preserve 
the magic symmetry. We consider the three different Z2 subgroup matrices (Eq. 4) of S3 as 
the perturbation matrices. Requiring symmetric mass matrices, the remaining S3 matrices 
are not considered because S< 123 ) and S< 132 ) are not symmetric and does not change 
the S3 invariance of the neutrino mass matrix. In Ref.j9] the breaking of S3 symmetry 
generated a non-zero #13 but the deviation from zero was small. Moreover, the perturbation 
parameter was considered to be real for simplicity leading to the absence of Dirac-type CP 
violation. Here, we generalize that work by considering a complex perturbation parameter 
which leads to Dirac-type CP violation in the model and, also, generates a relatively large #13 
consistent with the T2K results. We take the charged lepton mass matrix to be diagonal. If a 
horizontal symmetry exists it must, simultaneously, be a symmetry of the charged leptons as 
well as the neutrinos before the gauge symmetry breaking. After symmetry breaking when 
the fermions acquire mass, the charged lepton and the neutrino mass matrices should be 
constrained by different subgroups of the symmetry group in order to have non-zero mixing. 
We consider S3 to be the residual group in the neutrino sector. For the charged lepton 
sector, Z3 symmetry with the representation diag(l,co,uj 2 ) where u = e l2n ^ 3 , can be taken 
as the residual symmetry which yields non degenerate diagonal charged lepton mass matrix 
[TT] . Before the spontaneous symmetry breaking both the left-handed charged leptons and 
neutrinos share a common doublet so they must be governed by the same horizontal group 
which in this case comes out to be A(54) [12]. One may easily find that the gauge interactions 
of leptons remain invariant under the full group A(54). 



2 S3 invariant neutrino mass matrix 

The neutrino mass matrix (Eq. 6) can be diagonalised by any unitary matrix with a single 
trimaximal eigenvector, for example the TBM matrix Utbm- This S3 invariant neutrino mass 
matrix has eigenvalues a, a + 36 and a respectively. This is contrary to the experimental data 
since the two degenerate eigenvalues correspond to the mass eigenvalues nil and m 3 while 
experimentally m x and m 2 have smaller mass difference. A possible solution to this problem 
was provided by Jora et al. [8] by introducing a CP violating phase a between the complex 
vectors a and a + 36. The phase a is adjusted to ensure equal magnitude of a and a + 36, 
leading to degenerate mass spectrum for S3 invariant M v at the zeroth order. The complex 
plane is oriented in a manner so that the parameter 6 is completely imaginary, (Fig.(l)). 
The complex vector a lies in the fourth quadrant and is given by 

a = \a\e- ia/2 . (8) 

The complex parameters a and 6 in terms of the real free parameter x are 
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Figure 1: Angle a between two equal length vectors a and a + 3b. 



\a\ = x sec(a/2), 
2 

|6| = -xtan(a/2). 



(9) 



Democratic charged lepton mass matrices of the form bD and diagonal neutrino mass ma- 
trices of the form al have been discussed earlier [13] in the context of S3 symmetry. The 
possibility a = has, also, been discussed in the literature [H]. S3 flavor symmetry in the 
context of Type (I+II) see saw has been discussed in Ref.[15]. A generic feature of this 
approach is a democratic charged lepton mass matrix and a diagonal neutrino mass matrix 
at the zeroth order. Democratic mass matrix in the neutrino sector with terms breaking 5*3 
symmetry to \i — r symmetry has been discussed in [16] . In contrast, we consider the most 
general S3 invariant neutrino mass matrices at the zeroth order in the flavor basis where the 
charged lepton mass matrix is diagonal. 

3 Deviations from S3 symmetry 

In this section, we introduce explicit perturbations to the S3 invariant neutrino mass matrix 
M v and study their effects on neutrino masses and mixings. We introduce the perturbations 
in terms of a small, complex and dimensionless parameter \e % ^ . The perturbed neutrino 
mass matrix upto second order perturbation is given by 



where M v is invariant under S3 and S p and S q can be any of S^ 12 \ S*- 13 -' and S*- 23 - 1 such 
that p 7^ q and /1 is the real parameter with dimensions of mass and magnitude of the order 
of one in the units of the absolute mass scale. The matrices S^ 12 \ S^ 13 ^ and S^ 23 -* given in 
Eq. (4) belong to the same class of S3. In this scheme of perturbation, the S3 symmetry 
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is broken by its elements in the same representation. Even after introducing the second 
order perturbation, the neutrino mass matrix holds the magic property which leads to a 
trimaximal eigenvector for the perturbed neutrino mass matrix M' v . In the present case, the 
first order perturbation is taken to be Xe^S^ [9]. This leads to the following neutrino mass 
matrix after the first order perturbation with Xe l ^S^\ 



MP 
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a + b + pXe** b 
b a + b 

b b + n\e** 



b 

b + /iAe^ 
a + b 



J 



This mass matrix (Eq. 11) now leads uniquely to TBM mixing but now the degeneracy in 
the neutrino masses is broken. The complex neutrino masses are given by 



mi — (x + //A costp) + i{[iX sin?/' — xtan(|)), 

m 2 = (x + fiX cos iji) + i(fiX smip + xtan(|)), (12) 
m 3 = (x + fiX cos iji) — i(fiX sin-?/; + xtan(|)). 

To the first order in A, the neutrino masses are given by 



| mi | = x sec | + fiX cos(~ + tp) , 
\m2\ = x sec | + fiX cos(f — ip), 
\m?\ = x sec | — fiX cos(| + ip). 



(13) 



The second order perturbation leads to deviations from TBM and generates a non-zero #13. 
This can be done either by adding X 2 e 2l ^S^ or A 2 e 2 *^S'^ 13 ^ as the second order perturbation. 
The two choices only differ in their predictions for t9 2 3 and ip. Here, we take X 2 e 2% ^ as 
the second order perturbation. The perturbed neutrino mass matrix after second order 
perturbation becomes 
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It is clear that M' still satisfies the magic property and thus has a trimaximal eigenvector 
(^3' 73' V3^ T P^J' ^ e corres P onc iing mixing matrix can be obtained by the application of 
a general 1-3 rotation to the TBM matrix from the right UtbmUis where 
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A(54) 
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Figure 2: Symmetry breaking chain. 



Now we discuss the symmetry breaking chain, (Fig. (2)). After spontaneous symmetry 
breaking, the full group A (54) breaks to Z 3 in the charged lepton sector and 5*3 in the 
neutrino sector. Note that the S3 invariant neutrino mass matrix M v possesses accidental 
magic symmetry (Z' 2 ) i.e. M v is also invariant under Z' 2 
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The first order perturbation breaks 5*3 to the Klein's four-group Z 2 x Z 2 , where Z 2 is the 
H — r exchange symmetry. Although Klein's four-group is not a subgroup of S3 but the S3 
invariant M u has additional Z' 2 symmetry which allows for this breaking pattern. Z 2 x Z' 2 
is the symmetry group for TBM in the perturbed neutrino mass matrix M^. The second 
order perturbation breaks Z2 x Z 2 to Z 2 and, thus, the perturbed neutrino mass matrix M' v 
has a trimaximal eigenvector. 

Diagonalizing the perturbed neutrino mass matrix M' v using UtbmUi^, the neutrino masses 
to the second order in A are given by 
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The neutrino masses are made real positive by the phase matrix 



(18) 
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r 



a 



P 



| tan 1 i 



-xtan(^ 



-A/i sin(^)- 
\2~ 



-^-/xsin(2V>)- 



x+X^t cos(V')— 4-^ cos(2V>) 



1 f -If staii(f)+A/jsiii(?/>)+A 2 sin(2^) N 
° 2 Lcl11 V 3 ;+Amcos(^)+A 2 ) ucos(2^) /' 

l tan' 1 ^ " x tan (f )~ A/i sin W+V a* sin(2^) - 



(19) 



Thus, the full mixing matrix is given by 



U — UtbmUisP. 



(20) 



The solar and atmospheric mass-squared differences to second order in A are given by 



Am 2 2l = {m 2 2 - m\) « xA//(3A cos(2^) + 2(2 + A cos(V>)) sin(^) tan(f )), 
Aml 3 = (ml — ml) m xA/i(A cos(2^) + cos(^)(4 + 6A sin(^) tan(f ))). 



(21) 



In the earlier analysis [9], the mass-squared differences were independent of the phase a 
implying that it is a Majorana-type phase as lepton number conserving neutrino oscillations 
do not depend on this phase. This was the consequence of taking the perturbation parameter 
A to be real. In the present study, however, both the phases a and if) appear in the above 
mass-squared differences and the phase a can no longer be identified as a pure Majorana 
phase. If we take the phase if) to be zero, a decouples from the mass-squared differences 
and its Majorana character is restored. The neutrino mixing matrix after the second order 
perturbation is given by 
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where 
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The neutrino mixing angles upto the second order perturbation are given by 
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Parameter 


(+lcr,+2cr,+3cr) 

mean/ 1 o \ 


Am 2 21 [10- 5 eV 2 ] 


7 ctq(+0.20,+0.40,+0.60) 
' • i:)J (-0. 18 -0.35 -0.50) 


Am 2 31 [lO~ 3 eV 2 } 


o p- n (+0. 09, +0.18, +0.26) 
Z - OU (-0.16,-0.25,-0.36)' 
/ 9 4f ,(+0.08,+0.18,+0.27)\ 
v Z,4U (-0. 09, -0.17,-0. 27)) 


sin 2 9 12 


D oi o(+0.017,+0.038,+0.048) 
U - di/: (-0.015, -0.032,-0.042) 


sin 2 #23 


n r 9 (+0.06, +0.09, +0.12) 
U - OZ (-0.07,-0.11,-0.13)' 
(C\ C9(+0-06,+O.09,+O.12)x 
v U ' OZ (-0. 06, -0.10,-0. 13)/ 


sin 2 6>i3 


D fll o(+0.007,+0.015, +0.022) 
U - UiJ (-0.005,-0.009,-0.012)' 

l U - Uit) (-0.006,-0.011,-0.015) J 



Table 1: Current Neutrino oscillation parameters from global fits [19] . The upper (lower) 
row corresponds to Normal (Inverted) Spectrum, with Am 2 x > (Am^ < 0). 



The CP violation in neutrino oscillation experiments can be described through a rephasing 
invariant quantity, J C p p2] with J C p = I^{U , el U' 2 U'^U , * 1 ). The expression for J C p to 
second order in A is given by 

A 2 /icos(f)sec(f + V0sin(V0 fOQ , 
J cp Wx • (28) 

Since Jcp contains two different phases neither of them can be directly identified with the 
Dirac-type phase 5 in the standard PDG representation [18], rather some combination of 
these phases corresponds to 5. As pointed out earlier we have two choices for second order 
perturbation matrices. The above expressions were obtained for A 2 e 2l ^S'^ 12 ^ as the second 
order perturbation. If we use X 2 e 2l ^S^ as the second order perturbation, the predictions 
for all the parameters remain the same except for #23 and Jcp- The atmospheric mixing 
angle is now given by 

^, faM i_A£«a±M«^i+jw (29) 

while Jcp has the same value as above but with an opposite sign. The atmospheric mixing 
angle has the same deviation from maximal mixing as in the earlier case but in the opposite 
direction. 



4 Numerical results 

The experimental constraints on neutrino parameters at 1, 2 and 3<r [19] are given in Table 
1. The effective Majorana mass of the electron neutrino M ee which determines the rate of 
neutrinoless double beta (NDB) decay is given by [20] 

M ee = \ mi U' 2 + m 2 U[ 2 2 + m 3 U[l | . (30) 



S 




Figure 3: Correlation plots for the Normal Spectrum. 

Part of the Heidelberg-Moscow collaboration claimed a signal in NDB decay corresponding 
to M ee = (0.11 — 0.56)eV at 95% C. L. [21]. This claim was subsequently criticized in [22J. 
The results reported in [21] will be checked in the currently running and forthcoming NDB 
experiments. There are a large number of projects such as CUORICINO[23], CUORE [24J, 
GERDA [25], MAJORANA [26], SuperNEMO [27], EXO [28], GENIUS [29] which aim to 
achieve a sensitivity upto O.OleV for M ee . Cosmological observations put an upper bound 
on the sum of neutrino masses 

3 

s = E m - ( 31 ) 

1=1 

The WMAP data [30] limit X to be less than 0.67eV at 95% C.L.. For numerical analysis, 
we use the 3a ranges of neutrino oscillation parameters. For Normal Spectrum (NS), the 
parameters A and fi are constrained to lie in the range (0.012 — 0.135) and (0.1 — 2)eV, 
x is constrained to the range (0.034 — 0.234)el^ and the two phases a and ip have ranges 
(0.5° — 10.5°) and (92° — 120°) respectively. The effective Majorana mass M ee has the range 
(0.043 — 0.232)eU while the CP violating parameter Jcp lies in the range (0 — 0.07). For 
Inverted Spectrum (IS), the allowed ranges for the parameters A and \x are (0.011 — 0.115) 
and (0.1 — 2)eV while x can vary in the range (0.032 — 0.233)e\^ and the phases a and ip lie 
in the range (1.4° - 9.35°) and (57° - 86°) respectively. M ee lies between (0.065 - 0.235)eV 
and Jcp lies between ((—0.07) — 0). The deviations of the mixing angles for NS and IS are 
shown in Figs. (3a, 3b) and Figs. (4a, 4b) respectively. 
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Figure 4: Correlation plots for the Inverted Spectrum. 

5 Summary 

In the present work, we generalized the approach used in Ref.[9] for breaking £3 symmetry 
in the neutrino mass matrix by considering a small but complex perturbation parameter. 
Our starting point is a zeroth order mass matrix arranged to give three degenerate neutrino 
masses. This mass matrix can be diagonalized by e.g. the TBM mixing matrix. Mass 
splittings are introduced by adding perturbations to the zeroth order neutrino mass matrix. 
This also leads to a relatively large reactor mixing angle and CP violation in the neutrino 
mass matrix. The perturbation matrices are chosen to be the S3 group matrices in the 
three dimensional reducible representation belonging to the same class. Replacing the real 
perturbation parameter used in the earlier work [9] by a complex perturbation parameter 
has profound consequences for the phenomenology of neutrino masses and mixings. A real 
perturbation parameter produced very small reactor mixing angle and there was no CP 
violation in the model. To accommodate the recent T2K results indicating a relatively large 
reactor mixing angle it becomes necessary to consider a complex perturbation parameter 
which also generates CP violation. The neutrino masses are found to be quasidegenerate and 
a large effective neutrino mass is predicted which is testable in the ongoing and forthcoming 
experiments for neutrinoless double beta decay. 
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